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Abstract 

The NMR relaxation data on Sr 2 Cu0 3 [Phys. Rev. Lett. 76,4612 (1996)] are 
reexamined and compared with the analytic theory of the dynamic susceptibil- 
ity in the S=l/2 antiferromagnetic Heisenberg chain including multiplicative 
logarithmic corrections [Phys. Rev. Lett. 78, 539 (1997)]. Comparisons of 
the spin-lattice and the gaussian spin-echo decay rates (1/Ti and I/T^g) and 
their ratio all show good quantitative agreement. Our results demonstrate 
the importance of the logarithmic corrections in the analysis of experimen- 
tal data for quasi-lD systems and indicate that the dynamics of S^CuOs is 
well described by a S=l/2 one-dimensional Heisenberg model with a nearest 
neighbor exchange. 
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The spin- 1/2 one-dimensional antiferromagnetic Heisenberg model 

n = j y £s i -s i+1 (i) 

i 

continues to reveal novel many-body quantum effects!! in spite of its simplicity and extensive 
studies over many decades. In particular, a good understanding of the dynamic properties 
at finite temperature has been obtained only recently, using field theoriesiBBi and advanced 
numerical techniques.il Experimentally, there are generally two complementary methods to 
observe spin dynamics. The spin correlation function over a wide frequency (a;) and momen- 
tum (q) range can be measured by neutron scattering experiments.! The nuclear spin-lattice 
and the gaussian spin-echo decay rates (1/Ti and 1/T 2 g) measured by nuclear magnetic 
resonance (NMR) experiments! provide accurate information on g-averaged low frequency 
dynamic and static susceptibilities, respectively. Since it is the static and low frequency 
{Tiuj <C J) properties at low temperatures (T <C J) that can be treated most accurately by 
field theories, NMR is particularly suitable to test them. 

A remarkable consequence of the field theories is the quantum critical scaling behavior,i 
i.e, scaling of the dynamic susceptibility u;,T) in the variables c(q — ir)/T and Tiui/T 
(c = 7rJ/2 is the spinon velocity). This implies a divergence of the antiferromagnetic correla- 
tion length and the characteristic time scale of spin fluctuations as 1/T at low temperatures 
and reflects the critical nature of the ground state. The scaling determines the temperature 
dependences of the NMR relaxation rates as 1/T\ = const and 1/T 2 g oc l/v^T.! This behav- 
ior has indeed been observed in the experiments on Sr 2 Cu03,i which has the most ideal one 
dimensional character known to date.tMJ The theory also shows that the scaling is correct 
only approximately and there is a multiplicative correction with logarithmic temperature 
dependence.00 However, including this correction to the highest order actually degrades 
the agreement between theory and the data on Sr 2 Cu03.H 

Very recently, an improved analytic theory of the logarithmic corrections has been pro- 
posed by Starykh, Singh, and Sandvik,! generalizing the scaling ansatz and the conformal 
mapping to include the logarithmic factors. In this paper, we reexamine the NMR data on 
Sr 2 Cu03 reported in Ref. [8] and compare them with this theory. In the following, we first 
describe the theoretical results for the NMR relaxation rates both with and without the 
logarithmic corrections. We then reanalyze the spin-echo decay data to extract 1/T 2 g more 
accurately, taking the nuclear spin fluctuations due to spin- lattice relaxation into account, 
and also examine the appropriate value of J in Sr 2 CuC>3. From the comparison, we found 
that including the logarithmic corrections considerably improves the quantitative agreement. 

Let us start with the expression for 1/tM and l/T 2 Gc3t-J due to the magnetic hyperfine 
interaction 

H hi = ^ahaSja (a = a, b, and c) (2) 

between a nuclear spin at i site and an electron spin at j site, 

d Q ( a2,„\ , ,12/^1 Imx(g,^o) 

UJ 



T, 



2G 



1 






h 2 




P 


)'- 




8h 2 



f^{Al( q ) + Al( q) }^^> (3a) 



(3b) 



2 



Here x(q> u ) (x(q)) is the dynamic (static) susceptibility per spin in units of (g/is) 2 , u>o is 
the NMR frequency, and A a (q) = J2j A l £ ex P(.i<l r ij)- The magnetic field is assumed to be 
applied along the c direction. The prefactor in Eq. (3b) is chosen for nuclei with spin 3/2 
and is valid when the spin-echo decay is measured on one of three NMR lines split by the 
quadrupole interaction. The isotopic abundance p is 0.69 for 63 Cu nuclei. Since x(o) and 
lmx(q,w) are strongly peaked at q = it (the antiferromagnetic wave vector), A a (q) in the 
integral can be replaced by A a (ir). 

An analytic result for x(g, u) at finite temperatures was first obtained by Schultzjl using 
the bosonization technique to transform the Heisenberg model to a free boson Hamiltonian. 
The result satisfies the quantum critical scaling, 

x(q,oo) = Df 1 (q,u)/k B T J (4) 

where q = c(q—ir) / k B T, Co = hu/ksT, and D is an unknown constant determining the overall 
magnitude of x(g, u;). It immediately follows from Eqs. (3) and (4) that 1/T\ = const, and 
I/T2G oc 1/VT. In order to eliminate dependence on material parameters, we define the 
following normalized dimensionless NMR rates, which should be universal functions of T/J 
and can be compared directly with theories, 

(V^Onorm = ( TTTT ~, 7\ ( 5a ) 
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Sachdev used the result by Schulz to calculate the NMR rates and obtained! 

(l/71) noro = 2^, (VT/T^)^ = 1.1908D, (6) 

with their ratio being a universal number {T 2 G/Ti\/T) noTra = 1.680, which can be tested 
experimentally. The T-dependence of \/T\ and 1/T 2 q and the value of T 2 g/(VTTi) indeed 
agree quite well with the experimental data on Sr 2 Cu03.i 

However, the free boson theory neglects the marginally irrelevant operator in the original 
Heisenberg Hamiltonian, describing umklapp scattering processes, that leads to a multiplica- 
tive logarithmic correction to x(g, u). Sachdev has shown that to the leading order both 
(l/Ti) norm and (VT /T 2 g) norm acquire an identical multiplicative factor of ln 1 ^ 2 (T /T), where 
To is the high energy cutoff of the order of J.i However, the agreement with the experimen- 
tal data becomes worse if the factor /T) is included! Recently Starykh, Singh, and 
Sandvik have proposed a more elaborate theory of the logarithmic corrections! by adopting 
a simple ansatz generalizing the finite-size scaling and the conformal mapping to correlation 
functions with multiplicative logarithmic factors. They showed that the susceptibility takes 
the form, 

X (q, u)=D ln^To/T)^?, Co, A)/k B T, (7) 

with A = (1/4) {1 — 2in(T /r) }• m addition to the multiplicative factor, the logarithmic 
T-dependence of A also breaks the scaling. The equal-time real-space correlation function 
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obtained from this formula agrees very well with Quantum Monte Carlo (QMC) result for 
8 < J/T < 32 if Tq/ J ~ 4.5 and D « 0.075. The theoretical expressions for the NMR rates 
are! 

= ZD x 2 5 / 2 - 2A sin(27rA)/ 1 (A) In 1 / 2 (T /T)/tt 2 (8a) 

(VT/T 2G ) noTm = D x 2- 5 / 2+2A sin(27rA)r 2 (l - 2A)/ 2 (A) ln 1/2 (T /T)/7r 3/2 (8b) 

with 7 1 (A) = fo^dxj-^ and J 2 (A) =AJS°dx I^Syf- 

The above finite-T QMC estimate of the scale factor D compares well with recent T = 
numerical calculations (0.067890 and 0.0650), but certainly may contain some errors arising 
from non-asymptotic contributions in the temperature regime studied. The high-energy cut- 
off Tq has also be determined from fits of QMC data for the static staggered susceptibility and 
structure factor to their corresponding analytic expressions. ^ The results are T 0jX = 3.9 ±0.3 
from the susceptibility data and To iS = 5.1 ± 0.2 from the structure factor. Within the 
accuracy of the procedure these numbers can be considered being in good agreement with 
the value T = 4.5 quoted above. 

We now turn to the experimental results. NMR experiments on 63 Cu nuclei in Sr 2 Cu03 
have been performed by Takigawa et a/..i They have measured 1/Tj along the three crystal 
axes and I/T^g along the c direction. A static approximation was used to obtain 1/T 2 g 
from the spin-echo decay data, namely time dependence of I z due to spin-lattice relaxation 
process was neglected. Although this is a reasonable approximation in Sr 2 Cu03, for which 
(1/Ti) c is more than an order of magnitude smaller than (1/T 2 g) c , it is desirable to take 
the I z fluctuation effects into account to obtain more accurate values of 1/T 2 c. Recently, 
Cnrro et aB have derived an highly accurate analytic expression of spin-echo decay in the 
presence of the I z fluctuations, using the gaussian approximation proposed by Recchia et 
a/.J 2 ! We fitted the spin-echo decay data in Sr 2 Cu03 to this expression and extracted 1/T 2 g 
as defined in Eq. (3b). 

In order to convert the measured NMR rates to the normalized units defined in Eq. (5), 
we need the values of A a (Tc) and J. The values of A a (ir) / \/~J were determined from the width 
of the characteristic broad background of the NMR spectra due to a field-induced staggered 
magnetization near impurities,^ using the calculation of alternating local susceptibility near 
chain ends by Eggert and Affleck.! The accuracy of the these values depends on how closely 
the actual impurities (most likely the holes in the Zhang-Rice singlet statesi! due to excess 
oxygen) behave as chain ends. We consider it to be of the order of 10%. The anisotropy of 
A a can be determined much more accurately, with an uncertainty of a few percents. 

The exchange J was estimated to be 2200±200 K§ by fitting the temperature dependence 
of the uniform susceptibility to the theory by Eggert, Affleck and Takahashi.0 This value 
was used in the analysis of the previous NMR experiments.^ 2 ^ Suzuura et at, on the other 
hand, measured the optical absorption spectrum at T=32 K,@ which wasinterpreted to 
be due to simultaneous phonon and spin excitations. Lorenzana and Ederc^l analyzed the 
spectrum and obtained J=2850 K. Since this value is directly determined from the sharp 
peak of the spectrum located exactly at 7rJ/2 (the top of the des Cloizeaux- Pearson mode) 
except for a small shift corresponding to the optic phonon frequency (~ 0.08 eV), this latter 
value should be more reliable. 
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We speculate that the discrepancy may be resolved if one takes into account lattice 
degrees of freedom. A temperature dependence of J naturally results from the thermal 
lattice expansion. In fact, we recognize a slight systematic deviation between the x(T) data 
and the theoretical curve assuming a constant J (Fig. 4 in Ref. |ll|). Since the measured 
susceptibility includes the constant orbital and diamagnetic contributions, which are not 
known well, we can choose J = 2850 K at T = and allow J to change with temperature to 
reproduce the xCO data. We found that the x(T) data can be explained if J is reduced to 
2530 K at T = 800K, which is the highest temperature of the measurements. Since J oc t 4 
and t oc d~ 3 ' 5 , where t is the Cu to O transfer integral and d is the Cu-0 distance,^ thermal 
expansion of 0.8 % over 800 K is enough to account for this change of J. This seems quite 
plausible. 

A recent QMC study of the Heisenberg chain including dynamic (fully quantum me- 
chanical) phonons indicates that the fluctuations in J may furthermore lead to an apparent 
shift of J as obtained from x(T). A fit to the Heisenberg x(T), assuming a constant J, 
gives a result which is lower than the actual average spin-spin coupling in the spin-phonon 
model This effect may be the reason for the different values of J obtained from x(T) and 
the optical absorption. We use J = 2850 K in the following analysis but show also the 
results for J = 2200 K to indicate how sensitively the results depend on the value of J. 

Finally, we are at a position to compare the experimental results with the theory. Fig. 1 
shows the results for (l/Ti) nonn . Since the ratio A a (ir)/y/j is determined from experiments, 
(l/Ti) norm is independent of the choice of J. The solid line shows the result of Eq. (8a) 
with To/ J = 4.5 and D = 0.062, while the dotted line is the lnV 2 ( J/T) dependence with a 
magnitude chosen arbitrarily. It is clear that the theory by Starykh et al. shows much weaker 
T-dependence than ln 1 ^ 2 (J/T) and is closer to the experimental data. The value D = 0.062 
agrees with the value D = 0.075 determined from the comparison with the QMC results of 
the real-space correlation function! within the uncertainty due to errors in A a (jt)/y/~J and 
numerical determination of D. 

Fig. 2 shows the results for (vT/T 2 g) norm- The open (filled) circles are the experimental 
results for J = 2850 K (J = 2200 K) adjusted for the I z fluctuations. The crosses are the 
results for J = 2850 K in the static approximation. The correction for the I z fluctuations 
reduces the value of 1/T 2 q by 6 - 12 %. The solid and the dotted lines are the theoretical 
results of Eq. (8b) for D = 0.062 and D = 0.07, respectively (T / J = 4.5). We remark that 
only the first term in Eq. (3b), which we call the "scaling part", can be calculated reliably 
by the field theory, giving the results of Eq. (8b). The second term in Eq. (3b) is the square 
of the local susceptibility, which cannot be handled properly by continuum effective field 
theory.! It is smaller than the scaling part by a factor (T/ J) W(f /T). QMC calculations!! 
are useful to examine the relative magnitude of the second term. The filled and open triangles 
in Fig. 2 show the QMC results for the full and the scaling part, respectively. The QMC 
results for the scaling part agree very well with the analytic results with D = 0.07. The 
analytic results for D = 0.062 (the best value to fit the 1/Ti data) are in good agreement 
with the experimental results for J = 2850-fT (filled circles) if the negative contribution from 
the second term in Eq. (3b) (the difference between the two sets of QMC results) is taken 
into account. 

Fig. 3 shows the ratio of Eqs. (5a) and (5b). We used the 1/Ti data for H || a and H || b, 
since the hyperfine form factor A^,(q) + Ay(q) for H \\ c is relatively small at q = n. This 
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ratio provides a particularly stringent test for theories, since uncertainties in the values of 
D and A a (7r) cancels out and the only parameter that needs to be known is the exchange J. 
The results of the analytic theory are between the two sets of experimental data for different 
values of J. It predicts a weak logarithmic T-dependence, consistent with the experimental 
data, with an infinite slope at T = 0. The value at T = shown by the square in Fig. 3 is 
equal to the scaling result 1.680, which is temperature independent. Once again, including 
the logarithmic corrections improves the agreement with the experimental data. 

In summary, 1/71, V^2G, an d their ratio in Sr 2 Cu0 3 all show good quantitative agree- 
ment with the results of the analytic field theory that takes the logarithmic corrections into 
account. Our results demonstrate the importance of the logarithmic corrections in analyzing 
experimental data for quasi-lD systems, and that the dynamics in S^CuOs is well described 
by S=l/2 ID Heisenberg model with a nearest neighbor exchange. We suggest that remain- 
ing small discrepancies, such as the different estimates of J from the magnetic susceptibility 
and optical absorption experiments, may be due to spin-phonon interactions. 

We would like to thank S. Sachdev for interesting discussions and C. P. Slichter for 
communication on the analysis of spin-echo decay. This work was supported by the National 
Science Foundation under Grants DMR-89-20538 (AWS) and DMR-96-16574 (RRPS). 
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FIGURES 



FIG. 1. Results for (l/Ti) norm . Symbols are the experimental data for different field directions. 
The solid and dotted lines show the analytic theory Eq. (8a) with D = 0.062 and the ln 1//2 (J/T) 
dependence, respectively. 

FIG. 2. Results for (Vf/T 2G ) norm- The filled and open circles are the experimental data for 
different values of J corrected for the I z fluctuations. The crosses are the data for J = 2850 K 
without the corrections. The filled and open triangles show the QMC results for the full and the 
scaling (the first term in Eq. (3b)) parts, respectively. The results of the analytic theory Eq. (8b) 
for the scaling part is shown by the solid (D = 0.062) and the dotted (D = 0.07) lines. 

FIG. 3. Results for the ratio of (l/Ti) norm and (VT/T2G)noim- Symbols are the experimental 
data for different values of J and field directions. The result of the analytic theory (the ratio of 
Eqs. (8a) and (8b)) is shown by the line. Its value at T = (square) coincides the scaling result. 
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